Halo abundances and shear in void models 
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We study the non-linear gravitational collapse of dark matter into halos through numerical N-body 
simulations of Lemaitre-Tolman-Bondi void models. We extend the halo mass function formalism 
to these models in a consistent way. This extension not only compares well with the simulated 
data at all times and radii, but it also gives interesting clues about the impact of the background 
shear on the growth of perturbations. Our results give hints about the possibility of constraining 
the background shear via cluster number counts, which could then give rise to strong constraints on 
general inhomogeneous models, of any scale. 



PACS numbers: 98.80.-k 



I. INTRODUCTION 



Lemaitre-Tolman-Bondi (LTB) void models have been 
proposed as a viable alternative to dark energy. In these 
models one considers the possibility that we might live 
inside a large underdense region (a void) and that the 
apparent accelerated expansion of the Universe is only 
due to a misinterpretation of the observations in terms 
of a homogeneous background [51 [3H [33J in which the 
expansion rate is the same everywhere. The coincidence 
problem of ACDM (why now?) is substituted in these 
models by a violation of the Copernican Principle (why 
here?), since, in order to accommodate the observational 
constraints coming from the isotropy of the CMB and 
the matter distribution, the position of the observer is 
restricted to be very close (~ 1%) to the center of a 
highly spherical void. However, their ability to explain 
away many evidences for dark energy without any dark 
component has made LTB models a very attractive pos- 
sibility. 

It has been shown [H [9l [10] that a gigaparsec-sized 
void can reproduce reasonably well the distance-redshift 
relation deduced from current type la supernovae data. 
However, when these are combined with other cosmolog- 
ical probes, LTB models run into trouble. In particular, 
these models tend to predict a very high kSZ effect, due 
to the background contribution [TTJ [331 [35] ■ Also mea- 
surements of the local expansion rate combined with the 
full CMB power spectrum seem to be incompatible, the 
former being too low in LTB models [4, 23 , and more re- 
cently it has been shown that the latest BAO and SNe-Ia 
data show some tension too [37]. 

Nevertheless, research along these lines has been very 
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fruitful: these models have made us reconsider a non- 
standard approach to cosmology, using the machinery 
developed around them we have been able to consider 
observational effects due to the presence of large voids 
in a ACDM cosmology [TH1 (HI HH] , and we now under- 
stand much better the evolution of perturbations in an 
inhomogeneous background [7J [5U [35] . 

In a previous paper [5] we presented the first N-body 
simulations of LTB models, focusing on ensuring that 
the background evolution was correctly reproduced. In 
the present work we have performed the first study of 
halo statistics of one of the higher resolution simulations, 
and give predictions regarding the non-linear accretion of 
dark matter halos in void models. We will focus solely 
on the mass function and present a simple modification 
to the Press-Schechter theory allowing to accommodate 
and explain the effects introduced by using an LTB void 
model. 

The main interest of this approach is not just the char- 
acterization of halo abundances in LTB models, but a 
means to potentially distinguish whether the background 
space-time is FRW or not. Here dark matter halos are 
acting as probes of the growth of density perturbations 
and their mass function is extremely sensitive to the pres- 
ence of a finite background shear in large inhomogeneous 
voids of the LTB type. 

Moreover, we believe that smaller voids created due 
to the usual non-linear gravitational collapse associated 
with the cosmic web must also induce similar effects in 
the halo mass function, although at much smaller scales. 
If the approach used in the present work were applicable 
to these smaller voids, the reported contribution from 
the background shear could be included in the study of 
environmental effects on halo formation [141 122) . 
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II. THEORY 



The LTB metric 



The Lemaitre-Tolman-Bondi metric describes spaces 
with maximally symmetric (spherical) 2-dimensional sur- 
faces, and is given by 



ds z 



-dt z 



A' 2 (t,r) 
1 - k(r) 



dr 2 + A 2 (t,r)dn 2 , (1) 



for a matter source with negligible pressure and no 
anisotropic stress (X 7 ^ = — Pm(*> 1")Sq The function 
A(t,r) acts as an r-dependent scale factor. It is easy 
to see that in this framework the rates of expansion in 
the longitudinal (r) and transverse (#, <f>) directions are, 
in general, different (H T = A/A, H L = A 1 /A 1 ). With 
this setup, the Einstein equations can be written as an 
effective Fricdmann equation for a fixed r: 



H 2 (t,r)=H 2 (r) 



n(r) 



A 3 (t,r) 



+ (1 - fi(r)) 



A 2 (t,r) 



(2) 

where A (r) = A(t ,r) can be gauged to Ao(r) = r, 
Hq{t) = ffr(to,r) and O(r) is the ratio between the av- 
erage matter density inside a sphere of radius r and the 
critical density at that radius, and acts as an effective 
r-dependent matter parameter. 

The density profile of our simulated void follows the 
constraincd-GBH model jT0]. In it the free function fi(r) 
is parametrized by the central underdensity fl- ln , the void 
radius ro and the width of the transition void-background 
Ar/ro- The other free function Ho(r) is fixed by requir- 
ing a homogeneous Big Bang. This means that the void 
is a pure growing mode that disappears at very high red- 
shift [35]. We fix the underlying Friedmann- Robertson- 
Walker cosmology outside the void to be Einstein-de Sit- 
ter. For a more thorough discussion of the LTB metric 
and the constrained-GBH model, we refer the reader to 

0E5]. 



B. Linear perturbations and shear 

Although the equations describing perturbation the- 
ory in LTB models exist [7], we do not yet have a good 
understanding of their implications. The main prob- 
lem arises from the background not being homogeneous, 
and therefore the perturbations cannot be split into ir- 
reducible representations of 50(3) with decoupled equa- 
tions. In particular, the scalar modes couple to the vector 
and tensor modes through the background shear tensor. 
The background shear can be quantified as a normalized 
shear-to-expansion ratio 
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Ht — Ht 



2iJ 7 



Hr 



(3) 




FIG. 1: Normalized shear parameter e for the simulated LTB 
model. 



where E 2 = E^E 1 - 7 is the square of the background shear, 
and the expansion parameter of a congruence of comov- 
ing geodesies (see for further details). For voids of 
practical interests this shear parameter is usually small 
(in particular for the simulated model e < 0.06, see fig. 
[I]). Hence one would think that background shear effects 
can be neglected, in which case the equations for the den- 
sity perturbations reduce, at a fixed radius r, to those of 
an FRW universe with the corresponding effective cos- 
mological parameters at that r. We write this solution 



5 (t,r)<x D(n(r),A(t,r)/A (r)), 



(4) 



where 



D(n,a) = jWir(n,o) f - n d , a ' 

2 J [a' H(il, a') 



= a ■ 2-Fi 



(5) 



is the growth factor in an Open CDM universe with mat- 
ter parameter f2, with 2-Fi(a, b; c; z) the Gauss hypergeo- 
metric function. The full density perturbation equation 
has the form [3 



5 + 2H T 5 + [AH T + 6H^)S = O (S 2 , S 2 



(6) 



In the small shear limit, S — > 0, we propose the following 
Ansatz as an approximate solution 



S a (t,r) = 5 (t,r) (l + ae(t, r)Y 



(7) 



where the parameter a could, in principle, depend on 
(t, r) and possibly also on the cosmological parameters. 



C. The mass function in LTB models 

Most of the information about the non-linear accre- 
tion of dark matter halos is encoded in the mass func- 
tion n{M) dM: the comoving number density of halos 
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with mass M G (M, M + dM). The first theoretical de- 
scription of the mass function was developed by Press 
and Schechter .27] (PS hereon) and later re-derived and 
extended by Bond et al. in the so-called excursion set 
formalism [5]. Within this framework the abundance of 
halos can be predicted as the abundance of points in 
space in which the linear density contrast 5 smoothed 
over a scale corresponding to the mass M has crossed 
the spherical collapse threshold S c — 1.686. Although the 
PS prediction describes qualitatively well the mass func- 
tion, it fails to reproduce its details (overpredicting the 
density of low mass objects and underpredicting massive 
ones). Nevertheless it is a remarkable achievement that 
one can estimate the abundance of non-linear structures 
using only linear perturbation theory and the assump- 
tion that 6 is gaussian-distributed. The PS formula has 
been perfected using ellipsoidal collapse and empirical 
parametrizations [28] [31], so that n(M) can be calcu- 
lated to very good accuracy, often using one of the main 
results from this formalism: the mass function should 
be a universal (cosmology-independent) function of the 
variance of the linear density contrast field a(M, z) [T5] , 
Here we will use 



n(M, z) 



Pm 
M 



9(°) 



d In a 



dM 



(8) 



where a = er(M, z) and g(u) is given by |26j : 

, , abv b + 2cv 2 (l + ai> b ) . ,. 

9{(J) = (l + av")* &M ~ CV } ' 

(a, b, c) = (1.529, 0.704, 0.412), v = 5 c /a. 

We have also tried other parametrizations of the mass 
function j^Hl US] and checked that our results did not 
depend significantly on this choice. 

We follow the same rationale in order to calculate the 
mass function of halos at a given r and t in an LTB model: 
since the simulated void arises from a purely growing 
mode (i.e.: the Big Bang time is homogeneous), and per- 
turbations grow in a self-similar fashion, it is reasonable 
to assume that, in order to calculate the variance of 5m 
at (t, r), we should rescale the variance cr out (M, z) of the 
density perturbations outside the void, by a factor 



f(t,r) 



5 a (t,r) 



5 a (t, r — > oo) ' 



(9) 



where the density contrast is computed theoretically ac- 
cording to Eq. Q, and evaluated in Fig. [2j with and 
without the shear correction. Thus, our model for the 
mass function n(M, z, r) at a given radius r is Eq. Q 
with cr(M, z) substituted by <7 out (M, z) f(t, r). 

Note however, that our main results will be quoted in 
terms of the cumulative mass function within a sphere of 
radius r centered at the origin of the LTB patch: 



?i(> M, < r,z) 



4-7T r 3 



r' 2 dr' / dM'n{M',r',z), 
Jm 

(10) 



III. THE SIMULATION AND THE HALO 
CATALOG 

A more detailed description of the simulation we have 
used can be found in [2] (simulation V.). It has 960 3 par- 
ticles in a box of size L = 2400 Mpc hT 1 , which sets the 
mass resolution to m p = 4.2 x 10 12 M Q h . The simu- 
lated void has a size of tq = 1100 Mpc, a transition width 
of Ar/ rp = 0.3 and an underdensity of Q in = 0.25 (see 
The background cosmology is Einstein-de 
0.43. The small-scale pertur- 
1 and 
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Sitter (51 = 1) with h out 

bations are set using a power spectrum with n s 
erg = 0.9. It must be noted that an LTB void with these 
parameters is in fact ruled out, since the supernovae and 
baryon acoustic oscillation data seems to be mutually in 
conflict given a particular LTB profile [37]. We only use 
this simulation with a large background shear as a toy 
model to test our Ansatz about the mass function. The 
technique used to simulate LTB voids is based on setting 
the initial conditions (IC) appropriately by modifying the 
IC generator to take into account the large-scale pertur- 
bation induced by the presence of the void. This must 
be done at a high enough redshift so that the void can be 
regarded as a linear perturbation. A modified version of 
the 2LPT code [8] was used for this stage. Once the ICs 
are set, they are plugged into Gadget2 [SJJj, which we run 
in pure tree-mode. 

The halo catalog has been extracted using the AMIGA 
halo finder AHF [T3J [TTJ. It maps the particle content 
to an adaptively smoothed density field and locates the 
position of possible halos as local overdensities. Once the 
gravitationally bound particles around these have been 
extracted, the extent and mass of each halo is computed 



47T 

M(r vil ) = —r 3 vil . Pc A, 



(11) 



since this observable has better statistics. 



where we have used A = 200 as a collapse threshold [e.g. 

HBI. 



IV. RESULTS 

Once the halos have been identified, we can compute 
ri(> M, < r, z) merely by counting the number of ha- 
los with mass above M inside a sphere of radius r and 
dividing by the comoving volume of this sphere. The re- 
sults are shown in figure [2] It is easy to see (top panel) 
that, while approximating the growth of perturbations 
by its 0-th order in e (eq. El) yields a reasonably good 
fit at small and large radii (where the background shear 
vanishes), it fails to reproduce the halo abundances at 
intermediate radii. In the bottom panel we can see, how- 
ever, that adding a non-zero 1-st order correction (eq. [7]) 
solves this problem. Furthermore, we have found that 
this correction seems to be almost independent of r and 
t, with a ~ 2. 
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r (Mpc) 


&(« = 0) 


&(<* = 2) 


300 


0.093 


0.112 


600 


0.682 


0.104 


900 


8.321 


0.173 


1200 


2.394 


0.047 


1500 


1.474 


0.072 


1800 


0.478 


0.050 


2100 


0.177 


0.032 


2400 


0.095 


0.061 



TABLE I: Goodness of fit of our theoretical approach with 
and without a first order shear correction for different radii 
at z = 0. Similar results hold at all other redshifts (see table 

m. 



z 


&{* = 0) 


&(<* = 2) 





1.620 


0.072 


0.2 


1.163 


0.065 


0.35 


1.411 


0.042 


0.5 


0.667 


0.047 


1 


3.494 


0.824 



TABLE II: Goodness of fit of our theoretical approach with 
and without a 1-st order shear correction for different redshifts 
(summing over all radii). 



FIG. 2: Cumulative mass function for halos inside spheres of 
different radii compared with the theoretical prediction with- 
out the background shear correction (a = 0, top panel) and 
with a first order shear correction, where a — 2 V r, t (bottom 
panel). 



We have quantified the goodness of fit of our approach 
(with and without the shear correction term) using the 
measure: 



N 



'n(> M i} < r, z) ~ ?i 4 (r) N ~ 



n(> M i: < r, z) 



Here n(> M, < r, z) is given in section II C rii(r) is the 



cumulative mass function obtained from the simulation 
for a mass Mi within a sphere of radius r, and N is the 
number of mass bins. Figure [3] (top panel) shows the 
value of £jy for different choices of a at different r and 
z. A value of a ~ 2, found as the median of the best-fit 
values for all the calculated curves, gives a good fit in 
all cases with only a very mild dependence on r and t 
(shown on the bottom panel of Fig. |3|. The improve- 
ment due to the shear-correction term can also be seen 
in table [IJ in which we have calculated the goodness of fit 
with and without the shear correction for different radii 
at z = 0. Table [TT] shows the same result for different 
redshifts summing over all radii. This improvement is 
especially evident at intermediate radii, where e is larger 



and therefore its effects more important. 

It would be extremely interesting to investigate 
whether and how the value of a depends on the cos- 
mological parameters: if this parameter turned out to be 
independent of the void model, one should be able to pre- 
dict its value from some approximation in perturbation 
theory. However, this is work in progress and we defer 
the presentation of it to future work that will also make 
use of better resolved simulations. 



V. DISCUSSION & CONCLUSIONS 

We have extracted the halo content from an LTB N- 
body simulation and analyzed the halo abundances at 
different masses and radii. The main conclusions from 
this study are: 

• The theoretical description of the halo mass func- 
tion in FRW cosmologies can be fully extended to 
LTB void models by adding just one parameter (a) 
that accounts for the effect of the background shear 
on the evolution of matter density perturbations. 

• The value of this parameter (~ 2) seems to be con- 
stant in time and only mildly dependent on the 
position in the void. Whether this value depends 
weakly on the void model is still work in progress. 
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could potentially be used to constrain the amount of 
background shear, a crucial test for general inhomoge- 
neous cosmological models. A toy model that has been 
considered in the past in connection with LTB scenarios, 
and which could benefit from our analysis of background 
shear, is the swiss-cheese model [3l HH [20] . We leave for 
the future such investigation. It is also worthwhile ex- 
ploring whether our results apply to voids of astrophys- 
ical scales, of tens of Mpc, since they could have an ef- 
fect on the modelling of the environmental dependence of 
dark matter halo properties, as well as the backreaction 
of non-linear gravitational collapse on the background 
evolution. 




1200 1400 
r (Mpc) 



FIG. 3: Top panel: dependence of £ N (eq. 12 1 with the value 
of the shear correction term a at different r and t. Bottom 
panel: best fit value of a for different radii and redshifts, 
showing a mild evolution with r. 



If the shear correction turns out to be practically inde- 
pendent of the void model parameters, halo abundances 
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